A minimal model for kinetic arrest 
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To elucidate slow dynamics in glassy materials, we introduce the Figure-8 model in which TV hard 
blocks undergo Brownian motion around a circuit in the shape of a figure-8. This system undergoes kinetic 
arrest at a critical packing fraction cf> — cf>g < 1, and for <j) ~ 4>g long-time diffusion is controlled by rare, 
cooperative 'junction-crossing' particle rearrangements. We find that the average time between junction 
crossings tjc, and hence the structural relaxation time, does not simply scale with the configurational 
volume Q c of transition states, because tjc also depends on the time to complete a junction crossing. The 
importance of these results in understanding cage-breaking dynamics in glassy systems is discussed. 

PACS numbers: 64.70.Pf, 61.43.Fs, 82.70.Dd 
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Glass transitions occur in myriad systems that span a 
wide range of lengthscales [1] including atomic, poly- 
meric, and colloidal systems. When cooled or compressed 
sufficiently fast, glass-forming materials undergo a transi- 
tion from an ergodic liquid state to an amorphous solid-like 
glassy state |y,l3l01. Glassy dynamics is characterized by 
several common features [5]. For example, the viscosity 
and structural relaxation times diverge super-Arrheniusly 
near the glass transition, and the long-time self-diffusion 
constant becomes extremely small. Correspondingly, a 
plateau develops in the particle mean-square displacement 
(MSD). As the system approaches the glass transition, the 
plateau extends to longer and longer times fi H ISO, sig- 
naling kinetic arrest associated with the formation of cages 
of neighbors around each particle. Structural relaxation oc- 
curs through a series of rare cage-breaking events, in which 
particles in the nearest-neighbor and further shells move 
cooperatively so that caged particles can escape [4]. 

The last feature is emphasized in Fig. Q] where we show 
results from simulations of a 50%-50% bidisperse mixture 
of hard disks with diameter ratio 1 .4 in the supercooled liq- 
uid regime. In Fig. Q](a), we plot displacements of a focus 
particle as it moves between three cages a, (3, and 7. In 
Fig- (U(b)-(f), we monitor the focus particle and its neigh- 
bors as it breaks out of cage a and becomes trapped in 
/3. Two processes are required for cage-breaking to occur. 
First, an opening must appear in the shell of nearest neigh- 
bors (cf. panels (c) and (d)) and second, particles beyond 
the nearest-neighbor shell must make sufficient room to ac- 
commodate the escaping particle (cf . panels (d) and (e)). If 
the material beyond the nearest-neighbor shell behaves as 
a fluid, the second process occurs easily, and the first is the 
rate controlling step. However, at higher packing fractions 
(j) structural relaxation outside the nearest-neighbor shell 
relies on other cage rearrangements. As the system ap- 
proaches the glass transition, particles cannot escape from 
their cages, because the surrounding particles, trapped in 
their own cages, are unable to make room. Geometrically, 
this is similar to the mechanism responsible for rush-hour 
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FIG. 1: (a) Particle trajectory taken from molecular dynamics 
simulations of 64 bidisperse hard disks over a time period in 
which the focus particle explores three cages a, /3, and 7. The 
cage-entrance times are provided. The filled red circles corre- 
spond to the snapshots in (b)-(f), which are labeled by time. The 
focus particle is outlined in red and particles forming the a (/3) 
cage are outlined in blue (green). 

traffic jams, where cars cannot exit an intersection because 
there is not enough room in the street in front of them. This, 
in turn, prevents cars in the perpendicular direction from 
entering the intersection, which causes a cascade of delays, 
leading to a city-wide traffic jam. 

To elucidate geometrical mechanisms responsible for ex- 
tremely slow dynamics in glassy materials, we introduce 
a 'minimal' Figure-8 model that captures fundamental as- 
pects of caging and cooperative motion. As illustrated in 
Fig.lU in our model A^ hard blocks undergo single-file dif- 
fusion around a continuous course in the shape of a figure- 
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FIG. 2: Schematic of the Figure-8 model near 4> g . (a)-(c) de- 
pict a junction-crossing event. Initially half of the particles are 
in each lobe. Between (a) and (b) particle 5 crosses the junction. 
Between (b) and (c) particle 10 moves into the junction in the di- 
rection of the arrow. In (c) the center of particle 10 resides in the 
upper lobe, which completes the junction-crossing event. 

8. Kinetic arrest at large (ft occurs because particles mov- 
ing in one direction must vacate the junction to allow those 
moving in the perpendicular direction to pass through the 
junction. The Figure-8 model has several appealing fea- 
tures. First, it is one of the simplest continuum models that 
captures kinetic arrest (and the dependence of the structural 
relaxation time on (ft can thus be calculated analytically). 
Second, to mimic a glassy material, the Figure-8 model can 
be generalized to the 'Manhattan' model, which includes 
an arbitrary number of intersections and particles per lobe. 
Third, experimental realizations of the Figure-8 model can 
be performed, for examp le, by c onfining colloidal suspen- 
sions in narrow channels fl3l[l4ll . The Figure-8 model may 
provide insights into important unanswered questions con- 
cerning glassy systems: 1) What mechanisms give rise to 
cage -breaking events and how are they related to dynami- 
cal heterogeneities ijLLQI]? 2) Why does significant slow- 
ing down occur in dense particulate systems below random 
close packing? 3) What is the form of the divergence of the 
structural relaxation time near the glass transition? 

We now provide a detailed description of the Figure-8 
model. The upper and lower lobes of a channel of length 
L intersect at a square junction with unit dimensions. The 
particles (hard blocks with length I = 1) undergo single- 
file Brownian diffusion, implemented numerically using 
Monte-Carlo single-particle moves chosen from a Gaus- 
sian distribution. The total contour length L and aver- 
age gap A per particle satisfy the relation (ft = N/L = 
1/(1 + A). The particles are not allowed to turn at the 
intersection. Therefore, particles move through the junc- 
tion in one of two possible modes: (1) from north west to 
south east and vice versa or (2) from north east to south 
west and vice versa. To switch modes, particles moving in 
one mode must vacate the junction to allow particles in the 
other mode to enter the junction. We focus our analysis on 
small systems with 2 < N < 20 particles to imitate coop- 
erative cage rearrangements in finite-size local regions in 
glass-forming liquids (such as those depicted in Fig.Q]). 



The system dynamics is monitored by measuring the 

MSD of the blocks, S(r) = E^i(N( T ) ~ x t (0)] 2 )/N, 
where x,(r) is the position of the center of block i at time 
r. As seen from the numerical results shown in Fig. [3] 
(a), for sufficiently large (ft the MSD develops a plateau 
S = So, which signals the onset of slow dynamics. The 
length of the plateau increases as (ft approaches a critical 
value <fr g . However, for any (ft < <j) g , the MSD finally be- 
comes diffusive, with S(r) = 2Dr as r — > oo. This 
behavior mimics slow dynamics in glassy materials. The 
time td = S /2L> for the system to reach the long-time 
diffusive regime is plotted in Fig. [3] (b) as a function of (ft 
for several system sizes N. The results indicate that t d 
exhibits a power-low divergence 
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as indicated by the solid lines. 



To gain insight into the system dynamics near (fr g and 
explain the critical behavior CO, we consider sample tra- 
jectories depicted in Fig. |4] We see that most often par- 
ticles are evenly divided between the two lobes, and the 
junction is occupied by one or two particles in the same 
mode. Occasionally, the junction becomes unoccupied and 
the direction of motion changes. However, an analysis of 
the trajectories shows that not all such switches produce 
significant particle displacements. 

On rare occasions, a switch of the direction of motion re- 
sults in a significant shift of all particles. Such a junction- 
crossing rearrangement requires that a specific sequence 
of events unfolds as shown schematically in Fig. [2] and in 
actual trajectories in Fig. |4] First, a given particle com- 
pletely crosses the intersection and enters a lobe that al- 
ready has N/2 particles, thus creating a compressed lobe 
with N/2 + 1 particles. Next, another particle leaves the 
compressed lobe from the other end (which requires that 
the intersection is free) and enters the uncompressed lobe, 
so that the particles are again evenly distributed between 
the two lobes. Three such junction-crossing events are 
highlighted in Fig. |4](a), and a closeup of one of them is 
shown in Fig.|4](b). 

Due to geometrical constraints, junction crossings can 
only occur when N/2 + 1 < L/2 — 1, which yields 
(j) g = N/(N + 4). In the limit (ft — ► <f> g , the available 
space in the compressed lobe goes to zero, giving rise to 
the slow dynamics and the resulting plateau of the MSD 
depicted in Fig. [3] In the inset to Fig. [3] (a), we replot 
the MSD curves with r scaled by the average time be- 
tween junction-crossing rearrangements Tj C - The rescaled 
MSD curves for different (ft collapse at long times, which 
confirms that long-time diffusion is controlled by junction- 
crossing events. An analysis of the geometry of the system 
indicates that a junction crossing event produces an aver- 
age particle shift 5 = | + A — N^ 1 , and thus, we find that 
2D = 5 2 /t jc . 

Since the slow evolution for (ft — > (ft g results from an en- 
tropic bottleneck associated with creation of a compressed 
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FIG. 3: (a) MSD E(-r) normalized by A 2 for N = 10 at 
4> g — 4> = 0.1 (diamonds), 0.056 (squares), 0.032 (upward trian- 
gles), 0.018 (downward triangles), and 10 (circles). The time 
r is normalized by A 2 /2D S , where D s is the short-time diffusion 
constant. The dashed line is a fit to E(r) = 2Dt for large r at 
4> g — (f> — 0.032. The inset shows the same data versus time nor- 
malized by the average junction-crossing time tjc- (b) Diffusion 
time t d vs. <f> g - (p for N = 2, 4, 6, 8, 10, 12, 14 (from below). 
Fits to Eq. ((9} with a only weakly depending on TV (dashed lines) 
and asymptotic behavior (Q} (solid gray lines) are also shown. 

lobe, we expect that the scaling behavior of t d ~ Tjc with 
Ac/j can be obtained by calculating the corresponding vol- 
ume in configuration space. In equilibrium, the fraction of 
time t c /T the system spends in the compressed configura- 
tion can be expressed as 

i = if- <2) 

where Q cu and O are configuration integrals for the com- 
pressed (uncompressed) lobe and whole system respec- 
tively. fi c u and SI can be written in terms of the config- 
urational integral u(L Q , M) = (L - M) M /Ml for a ID 
gas of M unit-length hard rods confined within length L Q 
(Tonks gas [15]). The compressed and uncompressed lobes 
correspond to a Tonks gas of length L\ = L/2 — 1 with 
Nc, u = N/2 ± 1 particles. Accordingly we have 
^jv/2+i A^ 2 " 1 

^c = s^rL , TTT ; ^« = 7777^ 7T7 ; ( 3a > b ) 



(iV/2 + 1)! (iV/2-1)! 

where A e = \(L-N-A) and A„ = \(L-N) denote the 
free space in the compressed and uncompressed regions. 
The configurational integral for the whole system can be 
expressed as a combination of Tonks-gas results, 

n = 2co(L,N)- J2 w(.Li,JViM£i,JVa), (4) 

N t ,N 2 
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FIG. 4: (a) Trajectories Xi(r) scaled by A for each particle i 
over a time during which three junction-crossing events (high- 
lighted in orange) occur. The solid horizontal lines correspond 
to the center of the junction at x — (L/4, 3L/4, 5L/4, . . .), and 
the long-dashed horizontal lines are displaced by a particle size 
from them, (b) Close-up of particle trajectories during the first 
junction-crossing event. 

where the factor of 2 corresponds to the two directions 
of motion for particles in the intersection, and the sub- 
tracted sum prevents double counting configurations with 
an empty intersection. 

At (j) g the free space in the compressed lobe A c vanishes, 
which is the source of the kinetic arrest. The structural re- 
laxation time for <fi ~ (j) g is calculated using (0 to estimate 
the average time between junction crossings Tjc = T/n 
(where n is the number of crossing events). Assuming that 
on average during a single junction-crossing event the sys- 
tem spends time r c = t c /n with one compressed and un- 
compressed lobe, we obtain 
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(5) 



which links the structural relaxation time to the ratio of the 
configurational integrals (O and (O. Using ((3^), free space 
in the compressed lobe can be written in terms of A</j, 

VA _^/ 2 +! 



^~(0 s -0) yV/2+1 - (6) 

Thus, one might expect that t d y ~ O c ~ {<j> g - (j)) N / 2+1 , 
i.e. the inverse diffusion time scales with the number of 
transition configurations. (This assumption is frequently 
adopted in analyses of cooperative glassy dynamics [16]). 
However, our numerical results do not support this hypoth- 
esis, and instead we observe a weaker singularity ([]]). 

The anomalous behavior (fl]l stems from the fact that 
not only VL C but also r c vanishes at <p g . This can be 
demonstrated by noting that Brownian dynamics of a ID 
gas of hard rods can be mapped onto a ID system of 
point particles with interparticle distances equal to the 
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gaps Axi = Xi — Xi-x — 1 between rods in the origi- 
nal Tonks gas. Since a gas of point particles does not in- 
volve a characteristic lengthscale, the entire stochastic pro- 
cess (Axi(i), . . . , Ax M (t)) for systems of different pack- 
ing fractions (but corresponding initial conditions) can be 
scaled onto each other by introducing rescaled variables 

AxS) = X^AxiiX^t), (7) 

where A is an appropriate scaling factor (e.g. the average 
interparticle gap). In the above relations Ax\ is the gap 
between the first particle and the position at which the 
boundary condition is applied; the boundary condition on 
the other end of the domain is represented in terms of the 

variable ^- =1 Ax*. 

To apply scaling relation (0 to our system we recall that 
a junction-crossing event requires that the particles are di- 
vided between the compressed and uncompressed lobes. 
Particles in the compressed lobe evolve as a ID Tonks 
gas until an interaction occurs with a particle that initially 
resided in the uncompressed lobe. We note that the evo- 
lution of particles in the compressed lobe can be rescaled 
exactly even if a particle leaves this lobe and enters the 
junction, because both the particle positions and boundary 
conditions can be rescaled. This is important because a 
particle on the border of a junction enters and leaves the 
compressed lobe multiple times before a junction-crossing 
event is completed. 

The mapping (0 of the dynamics of the compressed lobe 
implies that the corresponding scaling will also hold for the 
average time t c that the system spends in the compressed 
configuration. Taking A = A c , we find that 

r c = aAf, (8) 

where a is a proportionality constant. Combining the 
above relation with (|5]l yields fl7ll 

_, aA- 2 n c tt u 



n 



(9) 



which, according to Eqs. (O and dU), agrees with the ob- 
served anomalous scaling behavior £[]). In Fig. [3] (b), we 
show that Eq. ([9} accurately represents the long-time diffu- 
sive dynamics not only in the scaling regime — ► <f) g , but 
also at moderate eft. 

To summarize, we introduced the Figure-8 model, which 
exhibits kinetic arrest when <p — > <j) g . We showed that for 
(j) ~ (f> g long-time diffusion in our system is controlled by 
rare, cooperative junction-crossing events, and we deter- 
mined the configuration-space volume £l c corresponding 
to the transition states associated with junction crossings. 
We also demonstrated that the inverse structural relaxation 



time r n does not scale with the volume S7 r as < 



, but 



the scaling also involves an additional singular factor asso- 
ciated with the accelerated evolution of compressed parti- 
cle configurations in the volume £l c . We predict that similar 
anomalous behavior may occur in glassy materials when a 
cage rearrangement requires compression of the material 
surrounding the cage. 



There are several possible extensions of the Figure-8 
model that may shed light on important features of the glass 
transition (such as dynamic heterogeneities and aging phe- 
nomena Illll[l2ll ). We expect that these phenomena can be 
characterized by expanding the approach used here B18I1 . 
One generalization of the Figure-8 model we are now pur- 
suing involves increasing the number of junctions j and 
determining how (j) g depends on j and N in a network of 
cross-linked channels. Similar to the Figure-8 model, the 
multi -junction system exhibits kinetic arrest at (j) g below 
close packing. An analysis of this class of models will thus 
shed light on the mechanisms that give rise to slow dynam- 
ics in glass-forming materials. 
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